Permutation-Symmetric Multicritical Points in Random Antiferromagnetic Spin 

Chains. 
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The low-energy properties of a system at a critical point may have additional symmetries not present 
in the microscopic Hamiltonian. This letter presents the theory of a class of multicritical points that 
provide an interesting example of this in the phase diagrams of random antiferromagnetic spin 
chains. One case provides an analytic theory of the quantum critical point in th e random spin-3/2 
chain, studied in recent work by Refael, Kehrein and Fisher (cond-mat/0111295). 
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Many of the interesting but poorly understood sys- 
tems of interest to current condensed-matter physics re- 
search are quantum many-body systems with both strong 
quenched randomness and strong interactions. One class 
of such systems where there are some experimental re- 
and significant theoretical progress has been pos- 
is antiferromagnetic Heisenberg spin chains. 
Much of the physics of these systems is captured by the 
model Hamiltonian 
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where the operator 5, represents a spin-S at site i of 
a linear chain. The nearest-neighbor exchanges J; are 
all positive, may be random, and may have an imposed 
"dimerization" 5: 



J i = J[l + S(-l) i ]e X p(Rr h ) , 



(2) 



where R measures the strength of the randomness and 
the rji are random numbers drawn from a distribution 
with mean zero, variance one, and all moments finite. 

This simple-looking Hamiltonian encodes a variety of 
low-energy behaviors depending upon the dimerization 
6, the randomness R, and the value of the spin S. For 
example, for S = 1, the undimerized chain (5 = 0) has 
a quantum critical point at some intermediate value of 
R that separates low-disorder Haldane and high-disorder 
Random Singlet (RS) ground states Hd]; ^ ms P om t is a 
multicritical point in the R-5 plane at which three phases 
meet ^ . In recent work, a related transition between low 
and high disorder RS states was also seen numerically in 
undimerized S — 3/2 chains B|. Here we examine these 
critical points in a larger context, showing that they are 
but two members of a countably-infinite class of random 
multicritical points. The low energy statistical properties 
of these special points exhibit the symmetry of the per- 
mutations of N identical objects, <Sjv, although for N > 2 
this is not a symmetry of the system's bare Hamiltonian. 

To proceed, we describe the phases of the spin chain 
in the valence-bond picture JTo| in which each spin-S 
is represented by the fully symmetrized multiplct of 2S 



spin- 1/2 's. The system's ground state has total spin 
zero (modulo end effects), so each such spin- 1/2 forms 
a singlet (a valence bond) with a spin-1/2 on a neigh- 
boring site. Thus, we can classify a ground state by 
how many such valence bonds are formed across the even 
links of the lattice: call this number a. Since each spin- 
s' must participate in 25 valence bonds, there must be 
(25 — a) valence bonds across the odd links. We will de- 
note this valence-bond solid (VBS) ground-state as being 
in the (cr,25 — a) phase, or, more compactly, the a phase. 
For spin-5, there are (2S+1) such phases, and various 
phase transitions between them (some phase diagrams 
are shown in Fig[j]). The VBS ansatz for the ground-state 
assumes that the valence bonds are all between nearest 
neighbors, which is not precisely correct even at R = 0. 
But there are indeed 25+1 topologically distinct pos- 
sibilities for the phases to which the real ground state 
can belong. As suggested by the valence-bond descrip- 
tion, these are distinguished by the properties of a chain 
end: For a chain in the a phase, if an even bond in an 
infinite chain is removed, the two resulting semi-infinite 
chains have ground states that contain free spin-(er/2)'s 
localized near their ends. 

The phase diagrams in the R-S plane are simple for 
zero or small R: At R — 0, all the phases (with a = 
0,1,... 25) can be accessed by sweeping S from -1 to +1, 
passing a succession of 25 critical points |T^]; for inte- 
ger 5, the (5,5) phase that occurs around S = is the 
familiar Haldane phase (Eg] , while for half- integer 5, the 
critical point between a = 5± (1/2) occurs at S — 0. At 
R = 0, the low-energy properties of critical points sepa- 
rating phases (a,2S~a) and (cr + 1,25 — a— 1) arise from 
residual spin-l/2's obtained by first forming er valence- 
bonds across the even links and 25 — a — 1 valence-bonds 
across the odd links. This leaves one unpaired spin-1/2 
per site, and these spins behave as a (critical) spin-1/2 
chain. Now, since the phases are gapped at R = 0, they 
survive for small R as well. Likewise, the critical points 
must extend to critical lines at R ^ 0, with the same 
low-energy properties as the random-exchange spin-1/2 
chain ||||: Along these lines, the chain is in the Spin-1/2 
Random Singlet (RSi/ 2 ) state. In this critical state, the 
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FIG. 1. Schematic phase diagram of S = 3/2 chains in the 
R-5 plane. Inset: Possible phases diagrams for 5 = 2 chains. 

residual spin- 1/2 at any given site "pairs" into a singlet 
with one at some other site not necessarily close to it, 
with the randomness determining the pairing; this pro- 
duces a glass of single valence-bonds with specific statis- 
tical properties at low energies ||. 

The behavior in the opposite limit of strong- 
randomness is also readily understood. Due to the broad 
distribution of exchanges at large R, the value of 5 is not 
crucial, and a spin-5 version (RSs) of the random singlet 
state obtains for 6 = 0; at low energies, such a state is a 
glass of 25-fold valence-bonds with the same statistical 
properties as the RS1/2 state. When 6 is non-zero and 
(say) positive, the valence-bonds in the RSs state will 
have their left end-point more often on an even site than 
on an odd site, and the resulting state is in the (25, 0) 
phase. Thus, turning on 6 in this regime drives the sys- 
tem into either the (25, 0) or (0, 25) phase, and the RSs 
state is the critical line separating these two phases. 

Given this picture of the phase diagram in the two 
limiting cases, one is immediately led to the interesting 
possibility that all the RS1/2 critical lines meet the RSs 
line at a single point at intermediate R and 6 = 0, pro- 
ducing a multicritical point at which all 25+1 distinct 
phases of the spin-5 chain meet. Indeed, a general the- 
ory of such multicritical points Vn at which N distinct 
phases meet forms the focus of the present Letter. 

We begin by addressing the question of existence: Re- 
cent work H shows that the 6 = transition from the 
gapless Haldane (1,1) phase to a RSi state in 5 = 1 
chains studied earlier [[>|J6|] is such a multicritical point 
(with N = 3), where all three phases of the system meet. 
In the recent RG study of the 5 = 3/2 case with 6 = 0, 
a single quantum phase transition between a RS1/2 state 
for small R and a RS3/2 state for large R was observed 
numerically M. Our discussion above shows that this 
transition is actually a multicritical point (with N = 4), 
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FIG. 2. A configuration of three adjacent domains. 

at which all four distinct phases of a 5 = 3/2 chain meet. 
For 5 = 2 or higher spin, the N = 25 + 1 multicritical 
point is not generically present when one only varies R 
and 6, to locate it requires tuning other parameters. For 
5 = 2, in particular, the possible topologies of the R-6 
phase diagram are shown in Fig 

Such multicritical points represent points at which the 
local ground state of the chain can be in any one of N 
possible phases. Thus, we develop a theory for the low- 
energy physics of these points in terms of the domain 
structure of the chain. We begin with a low energy pic- 
ture of the chain as being made up of a sequence of do- 
mains, each belonging to one of the N possible phases 
of a spin-5 chain (JV = 25 + 1). Neighbouring do- 
mains are separated by domain walls. These domain 
walls each carry spin: The magnitude of the spin on a 
wall is given by the number of unpaired spin-l/2's due 
to the difference of a across the wall. For two adjacent 
domains D\ and D2, with u\ and 02, as in Fig |2|, the 
spin on the domain wall separating them has magnitude 
5i2 = \o-\ - a 2 \/2. 

The low-energy properties of our chain are now con- 
trolled by the effective exchange couplings between neigh- 
boring domain- wall spins. In the absence of neighboring 
domain- walls, each domain- wall spin represents a zero- 
energy multiplet, with the spin localized near the wall. 
Neighboring domain-wall spins thus interact with an ef- 
fective exchange J that falls off rapidly with the domain 
length, and is consequently broadly distributed in mag- 
nitude (it can be of either sign). We allow each type of 
domain, a, to have its own probability distribution for the 
length of the domain and thus the exchange across the do- 
main. We thus have N probability distributions P a (j3\T) 
for the corresponding log-couplings (3 = ln(f2/| J|) > 0, 
where is the cutoff energy (the strongest exchange), 
and r = ln(f2o/f2) with fio a bare cutoff. 

The signs of the exchanges Ji in the domain picture are 
dictated by the domain sequence. Consider the configu- 
ration in Fig ||, assuming that J2 is the strongest of the 
three exchanges shown. At energy below \J 2 \ but above 
I Ji j and j J3 j , it should be possible to describe the system 
by replacing 5i2 and 523 with a single effective spin 5i3 
whose value is determined by the ground-state multiplet 
of the two-spin Hamiltonian J2S12 ■ S23. Consistency re- 
quires that this must be the same as eliminating D 2 and 
having a direct domain wall between D\ and D3 carrying 
spin 5i3 = |<7i — 03 1/2. For this to be true, J2 must be 
antifcrromagnetic (positive) if o\ — 02 and 03 — o<i are of 
the same sign, and ferromagnetic (negative) otherwise. 
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To proceed further, we need to specify the probabili- 
ties with which different domain sequences occur. We do 
this within a nearest-neighbor transfer matrix formalism. 
Thus we have a symmetric, purely off-diagonal N x N 
transfer matrix W aa i which gives the relative weights for 
the two types of domains a ^ a' to be present and adja- 
cent to each other. We normalize W to make its largest 
eigenvalue +1 (we denote the components of the corre- 
sponding normalized eigenvector as y/p^). This guaran- 
tees that the 'partition function' (sum over all possible 
domain configurations) Zl = Tr(W L ) for a sequence of 
L domains with periodic boundary conditions tends to 
unity as L — > oo. The unconditional probability for a 
given segment, say . . . p a PbPc ■ ■ •, to occur in the spa- 
tial sequence of domain-types is now given by modify- 
ing the expression for Z by introducing appropriate pro- 
jection operators n^ CT , = & aa i& ail at the corresponding 
places in the product of W's, yielding the expression 
Tr(. . . WW-WJl^WW-W . . .). Thus, the probability 
of the k th domain being type p is that of the k th 
domain being type (i and the (k + l) th being type v is 
y/PuPuWuy, etc. (for L — > oo). 

Given the broad probability distributions P a of the 
log-couplings (3 in our domain model, we can analyze 
the low energy properties using a strong- disorder RG ap- 
proach that eliminates, at each step, all excited states 
of the strongest-coupled pair of remaining domain-wall 
spins ||. For our domain- wall model, the RG action is 
rather simple, and in this strongly random limit does not 
generate any correlations between domains beyond those 
given by the nearest-neighbor transfer matrix W: Con- 
sider Fig and let J2 be the strongest bond. At each 
step, this RG "integrates out" the domain (in this case 
D 2 ) straddled by the strongest coupling: If o\ 7^ 03, 
D2 is eliminated, a direct domain-wall between D\ and 
D3 is formed, and the signs, but not the magnitudes, 
of Ji and J3 are altered (if necessary) to conform to 
the requirements of the sign-rule for this new configu- 
ration. If o\ = (T3, D2 is eliminated and D\ and D3 
are merged together into one domain. This merged do- 
main is straddled by a renormalized coupling of magni- 
tude IJ13I = IJ1J3/J2I and sign determined by our sign- 
rule. A little book-keeping yields the RG flow equations 
corresponding to this procedure: 
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V„(P a ® P a - P„) , 



dL 

df = - L[p 



P° + Y], 



(3) 



where V a0 = ^W a ^W^, 



P°. 



^(o|r), p a 

Pi - P2), Y 



p a = f J dfcdfaPMT^ifoirwp 

J2l> PvVw, P ■ P° = 2u P^Pp. an< ^ the sums run over the 
labels of the N domain types. Moreover, the flow of W 
also induces a change in p: 



^=-0 (P° 
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p-P°-Y). 
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We turn now to a fixed point analysis of these RG equa- 
tions. The multicritical point Vn is controlled by a fixed 
point with Sn statistical symmetry corresponding to 
freely interchanging between the N phases that meet at 
this point. This fixed point has W^ u = 1/(A^— 1) V/i^, 
P„(/?|r) = (iV-l)e-( iV - 1 )' 3 / r /r V/x, andp M = 1/JVV/i. 
Also, the number of domains decreases with the cutoff 
as L(T) = L(0)/r 1 /^ JV , with the exponent ^ N = 1/N. 
Thus, all N domain types are equally likely at this fixed 
point, with any two types of domains equally likely to be 
adjacent to each other. The fractions p s of the domain 
walls in the low-energy effective Hamiltonian that have 
spin-s follow simply from this. In the 5 = 3/2, N = 4 
case, we predict l/ip4 = 4, pi/ 2 = 1/2, pi = 1/3 and 
P3/2 = 1/6; the numerical estimates of [|| are in reason- 
able agreement with these results. [The low-temperature 
specific heat and susceptibility at the critical point are 
completely determined by ip an d Ps-] 

To get information on off-critical scaling properties, we 
need to analyze small perturbations about this Vn fixed 
point. Fortunately, the Sn symmetry imposes enough 
structure on the linearized flows to allow a full calcula- 
tion of all RG eigenvalues A which govern the T A growth 
or decay of the corresponding eigen-perturbations: There 
are only N — 1 relevant eigenvectors, all having eigenvalue 
A^r = (V4N + 1 - l)/2. Since one has to tune N - 1 
'knobs' in general to get N phases to all be 'degener- 
ate', this coincides with the minimal possible number of 
relevant directions at a AT- fold multicritical point — thus, 
this Sn fixed point governs all such strongly-random 
A-fold multicritical points. [In contrast, usual (non- 
random) multicritical points in Landau theory or in two- 
dimensions do not generically have Sn low-energy sym- 
metry. For example, at the critical fixed point of the 
two-dimensional three-state Potts model, there are four 
relevant modes, two of which produce flows to generic tri- 
critical points that do not have the £3 symmetry of the 
Potts model.] The relevant eigenvectors can be chosen to 
correspond to perturbations which make only one of the 
N phases fall out of favour, thus reducing the symmetry 
from Sn to Sn—i- In addition we have one irrelevant 
eigenvector with eigenvalue —1 (representing an additive 
shift in r), for N > 2 there are N —1 irrelevant eigenvec- 
tors with eigenvalue —(y/AN + 1 + l)/2, and for N > 3 
there are N(N — 3)/2 irrelevant eigenvectors with eigen- 
value —N. [We also expect || other 'infinitely' irrelevant 
perturbations, i.e. decaying exponentially with T — these 
are not considered here.] 

In the 5 = 3/2, A = 4 case, we thus predict a rel- 
evant eigenvalue of A J = (vTf — l)/2, and a correla- 
tion length exponent v = 1/A^V4 — 2.56 (note that the 
numerical estimate of Q differs significantly from this 
prediction, probably due to slow transients or finite-size 
effects). Deviations from Va in the R-5 plane contain 
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linear combinations of the three relevant perturbations. 
The fact that the RG eigenvalues are all the same means 
that the phase boundaries come in linearly at V4. The 
slope of two of the four phase boundaries is fixed by not- 
ing that any phase boundary between phases related by 
the <5 — s- — <5 symmetry of the problem (corresponding to 
an interchange of even and odd sites) must lie on the R 
axis. The same symmetry forces the the remaining two 
to be reflections of one another about the R axis. To 
fix their slopes, consider one of them, say the (2,1) to 
(3,0) phase boundary: These two phases are degenerate 
here, but the other two have higher energy. A positive <5 
lowers the energy of both phases (since they have more 
singlets on the even bonds than on the odd bonds), with 
the energy of the (3,0) phase lowered more than that of 
the (2,1) phase. Decreasing R, on the other hand, lowers 
the energy of the (2,1) phase relative to the (3,0) phase. 
Thus, these two phases will remain degenerate if R is de- 
creased with increasing 5, and this phase boundary must 
leave "P4 with a negative slope (see Fig. |l|) . 

The above RG equations also admit fixed points of 
lower symmetry Sm (with M < N) at which domains of 
M phases each occur with equal probability, and other 
domain-types do not occur at low energy. These fixed 
points govern locii of multicritical points Vm at which 
M of the phases meet. Two examples of such locii are 
the RSx/2 and RS 3 / 2 lines in the R-S phase diagram of 
5 = 3/2 chains (these have M = 2). Other examples 
(with M = 3) include the points V3 in the R-8 phase di- 
agram (see Fig [l]) of 5 = 2 chains — these have the same 
exponents as the point V3 of 5 = 1 chains. In addition 
to these lower-order multicritical points, the RG equa- 
tions also admit 'Griffiths' fixed points, which describe 
the continuously varying power-law singularities within 
the individual phases M. 

Turning to 5 > 3/2, it is now clear that the generic 
phase diagram in the R-S plane for 5 — 2 will look like 
one of the two insets shown in Fig [j], with the putative 1S5 
symmetric point "P5 splitting into lower-order multicriti- 
cal points as shown. All five phases of the 5 = 2 chain 
will only meet at V§ upon fine-tuning some additional 
parameter in the model (such as nearest-neighbor inter- 
actions more general than simply the exchange 5j ■ St+i). 
Similar considerations also rule out the generic occur- 
rence of such maximally symmetric multicritical points 
in the R-S plane for all 5 > 2. 

Finally, we note that the basic structure of the domain- 
wall model used here can also be motivated from a more 
microscopic argument with ([l]) as the starting point. 
Consider treating (|l|) for arbitrary 5 with a generaliza- 
tion of the approximate extended Ma-Dasgupta-Hu RG 
approach H^). [For general 5, the procedure eliminates 
all excited states of the most strongly coupled pair of 
spins if this coupling is ferromagnetic, while taking care 
to eliminate only the highest excited state if it is anti- 



ferromagnetic; the RG rules for signs and magnitudes of 
couplings are as in Ref ||[|.] Assign the formal domain 
label a = (a = 25) to every even (odd) bond of the 
unrenormalized Hamiltonian (]]]). Eqn (|l]) with this la- 
beling is consistent with the rules for domain-wall spins 
and signs of couplings in our domain- wall model; we can 
therefore formally think of each J of (|l]) as straddling a 
domain of the corresponding type. Now, if the renormal- 
ized Hamiltonian and choice of labels at a given stage of 
the RG is consistent with this domain interpretation, it 
is possible to relabel the couplings after each RG step 
to preserve this property: Let J2S12 ■ S23 in Fig || be 
the term with the largest gap between lowest and high- 
est energy states. If J2 > 0, and neither S12 nor S23 are 
spin-1/2, change a 2 to 02 + (cr\ — a^VIci — 0-2! after the 
next RG step (which reduces both S12 and S23 by 1/2). 
If J2 > and S12 — S23 — 1/2, attach the common label 
of J\ and J3 to the new coupling J13 that reaches across 
them after the next step (which puts them into a singlet 
state). In all other cases, there is no need to relabel any of 
the couplings that remain after the next step. Although 
this RG procedure is approximate, it is expected to be 
qualitatively accurate for low-energy properties ||||||, 
and the formal device above thus provides an alternative 
route to our domain model. 
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